Oscillations and interactions of dark and dark-bright solitons in Bose-Einstein 

condensates 
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Solitons are among the most distinguishing fundamental excitations in a wide range of non-linear 
systems such as water in narrow channels, high speed optical communication, molecular biology 
and astrophysics. Stabilized by a balance between spreading and focusing, solitons are wavepackets, 
which share some exceptional generic features like form-stability and particle- like properties. Ultra- 
cold quantum gases represent very pure and well-controlled non-linear systems, therefore offering 
unique possibilities to study soliton dynamics. Here we report on the first observation of long-lived 
dark and dark-bright solitons with lifetimes of up to several seconds as well as their dynamics in 
highly stable optically trapped 87 Rb Bose-Einstein condensates. In particular, our detailed studies of 
dark and dark-bright soliton oscillations reveal the particle-like nature of these collective excitations 
for the first time. In addition, we discuss the collision between these two types of solitary excitations 
in Bose-Einstein condensates. 

PACS numbers: 03.75.-b,03.75.Lm,05.45.Yv, 



The dynamics of non-linear systems plays an essential 
role in nature, ranging from strong non-linear interac- 
tions of elementary particles to non-linear wave phenom- 
ena in oceanography and meteorology. A special class of 
non-linear phenomena are solitons with interesting parti- 
cle and wave-like behaviour. Reaching back to the obser- 
vation of " waves of translation" in a narrow water chan- 
nel by Scott-Russell in 1834 [1 , solitons are nowadays 
recognized to appear in various systems as different as 
astrophysics, molecular biology and non-linear optics [2]. 
They are characterized as localized solitary wavepack- 
ets that maintain their shape and amplitude caused by 
a self-stabilization against dispersion via a non-linear in- 
teraction. While an early theoretical explanation of this 
non-dispersive wave phenomenon was given by Korteweg 
and de Vries in the late 19th century it was not before 
1965 that numerical simulations of Zabusky and Kruskal 
theoretically proved that these solitary waves preserve 
their identity in collisions [3j H]. This revelation led to 
the term "soliton" for this type of collective excitation. 
Nowadays solitons are a very active field of research in 
many areas of science. In the field of non-linear optics 
they attract enormous attention due to applications in 
fast data transfer. 

Bose-Einstein condensates (BEC) of weakly interacting 
atoms offer fascinating possibilities for the study of non- 
linear phenomena, as they are very pure samples of ultra- 
cold gases building up an effective macroscopic wave 
function of up to mm size. Non-linear effects like col- 
lective excitations [5 , four- wave mixing [6] and vortices 
[7J |H] have been studied, to name only a few examples. 
The existence and some fundamental properties of soli- 
tons have been deduced from few experiments employ- 
ing ultra-cold quantum gases. Bright solitons, charac- 
terized as non-spreading matter-wave packets, have been 
observed in BEC with attractive interaction [9j [10] [11] 



where they represent the ground state of the system. In a 
repulsively interacting condensate confined in a periodic 
potential, bright gap solitons have been realized p~2j [13] 
by modelling a suitable anomalous dispersion. 
BEC with repulsive interactions allow for dark soliton 
solutions, characterized by a notch in the density distri- 
bution. In contrast to bright solitons, dark solitons are 
truly excited states with energies greater than the under- 
lying BEC ground state. Dark solitons have been gener- 
ated in few pioneering experiments [14] [15] [16] boosting 
an immense theoretical interest in these non-linear struc- 
tures in BEC. Dynamical [TS1IIZ1IIH1IIS1I2D1EI1I221I231 
and thermodynamical [24] [25] [26] instabilities as well as 
collisional properties [27, 28] have been analyzed theoret- 
ically in great detail. Moreover, the existence of bright 
solitons stabilized by the presence of dark solitons in an- 
other quantum state in repulsive BEC has been proposed 
[22] and confirmed in a proof-of-principle experiment [16] . 
The occurence of undamped oscillations of solitons in ax- 
ially harmonic traps comparable to a single particle with 
a large negative mass has been predicted by mean field 
theory [22l \2E\ and supported by numerical simulations 
[24] for dark solitons and appears as one of the paradigms 
of soliton physics in BEC as it clearly demonstrates the 
particle character of the soliton. Limited by very short 
lifetimes in previous experiments however, oscillations of 
dark solitons have not been observed yet. 
In this paper, we show that very long-lived dark solitons 
can be generated in BEC, facilitating detailed studies of 
soliton oscillations as well as the creation of dark-bright 
solitons by filling the dark soliton with atoms in another 
hyperfine state. In addition the reflection of a dark soli- 
ton off a filled soliton has been observed and is presented. 
In the following we give a brief introduction to the de- 
scription of solitons in BEC and characteristic properties 
of soliton dynamics. Close to absolute zero temperature, 



where thermal fluctuations can be neglected to first or- 
der, the condensate is well described within the frame- 
work of the non-linear Gross-Pitaevskii equation (GPE) 
[29] , which is known to support soliton solutions as it 
is closely related to the cubic non-linear equation vastly 
employed in non-linear optics to describe solitary wave 
propagation in optical fibers [2] 

H 2 

ihil>(z,t)= - 2^" CM) 

+ \v ext {z)+g\^{z,t)\ 2 ^{z,t). (1) 

Here ip denotes the condensate wavefunction, g = 2hu±a 
is a measure of the non-linear atomic interaction in quasi- 
1D geometries with strong transverse confinement, g is 
determined by the s-wave scattering length a in quasi- ID 
and the transverse trapping frequency lj± . V ext describes 
a confining external potential. A dark soliton solution to 
the GPE of a homogeneous BEC describing a density 
notch at position q propagating along z with a velocity q 
can be written as [30] 
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where the speed of soun d in a quasi- ID condensate is 
given by c s = ^nog/2m. no is the peak density of the 
condensate. The inverse size of the soliton k, is deter- 
mined by the healing length £ = h/mc s and the soliton 
speed q through n = £ _1 x yl — (q/c s ) 2 . Note that 
in a quasi- ID condensate, the averaging of the density 
over the radial degrees of freedom effectively changes 
the speed of sound to c s = c s /V2 and the healing length 
to £ = >/2£ as compared to their 3D values c s and £ 
respectively. The phase and density distributions of 
a dark soliton are schematically shown in Fig. lb. As 
depicted, the phase only shows significant changes in the 
vicinity of the nodal plane of the soliton and is constant 
otherwise. Crossing the nodal plane of the soliton, the 
wavefunction accumulates a specific phase slip between 
and 7r depending on the depth and speed of the soliton 
related by n s /no = 1 — (q/c s ) 2 = sin 2 (0/2) where 
n s denotes the missing density at the position of the 
soliton. A phase jump of A0 = tt corresponds to a fully 
modulated soliton with zero velocity representing the 
only time independent soliton solution of the GPE. As 
the phase difference diminishes the soliton velocity grows 
while it gets shallower and wider ultimately vanishing at 
the speed of sound. The preceeding considerations have 
led to the idea and experimental realization of creating 
dark solitons in BEC by optically imprinting a phase 
gradient of approximately tt over a spatial region not 
larger than the healing length p~4l ISJ [31] . 
While solitons in quasi one dimensional BEC are dynam- 
ically stable p~8j [25] , experiments p~5j [16] and intensive 
theoretical studies suggest that in less constrained 



geometries, where the condition 7 = nog/hujj_ <C 1 is not 
perfectly met, the growth of dynamically unstable modes 
will lead to a transfer of kinetic energy of the soliton 
to radial excitation modes of the condensate, mediated 
by the atomic interaction and resulting in a bending of 
the soliton plane [15 , which may ultimately decay into 
vortex pairs as reported in Ref. [16]. Since the energy 
of a dark soliton is always greater than the energy of 
the ground state condensate it is thermodynamically 
unstable in any case and shows fast decay even at 
reasonable low temperatures [24j [26] [32] as observed 
in experiment [HI [15]. The dissipation accelerates the 
soliton according to its negative kinetic energy until it 
vanishes and smoothly transforms to the BEC ground 
state as it approaches the speed of sound [25]. This 
interesting aspect can be interpreted as an accelerating 
instability [22] and implies that a negative mass can 
be assigned to a dark soliton. Lifetimes on the order 
of 10 ms have been reported preventing the observation 
of more complex soliton physics like oscillations or 
collisions. 

We have developed a reliable, robust method to produce 
elongated 87 Rb-BEC at extremely low temperatures 
in an optical trapping potential overcoming former 
technical limitations. We produce a BEC composed of 
5 x 10 4 87 Rb atoms in the 5 2 Si/ 2 , F = 1, m F = -1 
state in an optical dipole trap with trapping fre- 
quencies uj z = 2tt x 5.9 Hz, uo Y f v = 2tt x 85 Hz and 
^hor _ 2^ x 133 Hz with no discernable thermal frac- 
tion. Trap frequencies have been cross-checked by the 
measurement of various collective oscillations. Typical 
atomic peak densities are 5.8 x 10 13 cm -3 implying a 
speed of sound of c s = 1.0 mm/s. The chemical potential 
is on the order of 20 nK. 

Ultra-stable laboratory conditions ensure an exceptional 
reproducibility enabling us to record time series of 
soliton dynamics with unprecedented precision. The low 
trap depth guarantees a slight but constant evaporative 
cooling, so that no heating can be detected for time 
scales as large as the lifetime of the condensate, which 
is greater than 10 s. Solitons are produced by optically 
imprinting a phase gradient as shown in Fig. la: a part 
of the condensate is exposed to the dipole potential 
Z7di p of a laser beam detuned by some tens of GHz from 
atomic resonance. We image an optical mask pattern 
onto the BEC with diffraction-limited optical resolution 
of better than 2 /am. This results in a phase evolution of 
the masked relative to the unmasked part of the conden- 
sate of Acj) = Udipt/(ih). The pattern is generated by a 
spatial light modulator (SLM) with an effective pixel size 
of 0.8 /im allowing for almost arbitrary optical potentials 
[33] , To imprint a phase slip of order n, we choose a pulse 
time t n = 40 /is, much smaller than the correlation time 
fcorr = £/c s = 700 /is for our experimental parameters to 
avoid a simultaneous disturbance of the atomic density. 
This phase gradient leads to a local superfluid velocity 
of the condensate according to vsf — h/md z <j) which 
can also be interpreted as a local potential gradient 



transferring momentum to the BEC, thus assisting the 
formation of a density minimum [14]. Since a dark 
soliton can be regarded as a hole rather than a particle 
it moves in the direction opposite to the superfluid flow 
of the condensate. The appropriate equation of motion 
for small velocities q reads M s q(t) = —l/2d z V(z) where 
M s is the negative mass of the soliton. This implies a 
soliton oscillation frequency of uj = co z /\^2 for harmonic 
confinement [22j [25] . 

Fig. 2a shows the time evolution of a dark soliton 
created by the aforementioned phase imprinting method. 
Absorption images were taken after a time-of-flight of 
11.5 ms to allow the condensate and soliton to expand 
since the soliton size l s ~ £ ~ 0.8 /im in the trap is 
beyond optical resolution. The soliton clearly propagates 
axially along the condensate with an initial veloctiy of 
q = 0.56 mm/s = 0.56 c s indicating a relative soliton 
depth of 0.68. We were able to detect nearly pure dark 
solitons after times as long as 2.8 s in single experimental 
realizations (see Fig. le), surpassing lifetimes of dark 
solitons in any former experimental realization by more 
than a factor of 200. Fluctuations in the soliton position 
due to small preparation errors however prevent the 
observation of soliton oscillations for evolution times 
Tevoi ^> 250 ms. The extraordinary long lifetimes 
facilitate the first observation of an oscillation of a 
dark soliton in a trapped BEC. An oscillation period of 
ft = 2tt x (3.8 ± 0.1) Hz has been recorded and could 
be followed for more than one period. In contrast to 
the theoretical prediction of u z /\/2 = 2tt x 4.2 Hz for a 
dark soliton in an axially harmonic trap the observed 
period indicates that the harmonic approximation for 
the trapping potential starts to break down. Caused 
by the shallowness of our dipole trap the atoms rather 
experience a Gaussian potential, which is less steep 
than harmonic leading to a larger amplitude-dependent 
oscillation period for the soliton. We have calculated the 
oscillation frequency in a Gaussian potential created by 
a laser beam with a waist of 125 /im with the observed 
soliton amplitude of Z s = 33 /im and find an oscillation 
frequency of 4.0 Hz. This is in good agreement with 
our experimental data. Furthermore, the observed 
amplitude allows for a consistency check of the soliton 
depth. At the turning point of the soliton motion Z s the 
constant soliton depth equals the Thomas-Fermi density 
utf(Zs) of the condensate and interrupts the superfluid 
flow of atoms. At this point the soliton starts to move in 
the opposite direction. Given the measured initial speed 
of the soliton and the observed density distribution of 
the condensate Z s can be calculated to be 36 /im and is 
in very good agreement with the measured value. 
Another feature extracted from Fig. 2a is a density wave 
that travels in the opposite direction at a velocity equal 
to the speed of sound. The occurrance of such a density 
wave has been investigated theoretically and experi- 
mentally [14] and has been attributed to the method 
of soliton generation via phase imprinting while leaving 
the instantaneous density distribution unchanged. The 



density waves die out after approximately 50 ms leaving 
a flat BEC with only one soliton excitation. 
Calculating the dimensionality parameter 

7 := nog/hco± = 3.7 and comparing this value to 
the critical ratio j c given by Muryshev et al. [26] we 
find our soliton to be right on the edge of the region 
of dynamical stability. This is confirmed regarding the 
observed soliton lifetimes. 

We have performed numerical simulations of the GPE 
showing that the phase imprinting method cannot create 
single perfect dark solitons but always creates density 
waves that carry away part of the imprinted phase 
gradient. Moreover the occurence of a second small 
soliton can be extracted from the simulations as shown 
in Fig. 2b. 

The crucial feature to the observed long lifetimes of dark 
solitons seems to be the very low temperature of our 
samples. The critical temperature for Bose-Einstein con- 
densation for our experimental parameters is (67 ±5) nK. 
Estimating that a thermal fraction of at least 10% could 
have been detected in absorbtion imaging -which was 
not the case- an upper limit for the temperature of 
T < 0.5 T c = 30 nK can be given which is on the order 
of the chemical potential /i. We assume a significantly 
smaller temperature, since temperatures of T « 0.2 T c 
would already considerably limit the solitons lifetime 
[24] . which has not been observed in our experiment. 
Despite the interesting physics that can be investigated 
using dark solitons, so called dark-bright solitons 
appearing in multi-component BEC show even more 
fascinating physical properties, such as enhanced dy- 
namical stability and the possibility of bright-component 
particle exchange in soliton collisions. A dark-bright 
soliton is basically a dark soliton filled with atoms of a 
different species or in another internal state of an atomic 
matter wave. A dark-bright soliton can be generated 
in a 8T Rb BEC by imprinting a dark soliton in state 
\F = l,ra.p = 0) and filling the density dip with atoms 
in state \F = 2 1 rriF = 0) [34] [35] leading to the density 
distribution depicted in Fig. 3b. While dark solitons 
are unstable to transverse excitations with wavelengths 
greater than their extension l s « £ large dark-bright 
solitons are expected to overcome this restriction since 
their size can be much larger than £ when the number of 
atoms in the other hyperfine state becomes very large. 
Dark-bright solitons should therefore be robust in traps 
geometries, which are not truly ID. We have employed 
a method to simultaneously imprint the phase gradient 
and transfer atoms to the other hyperfine state via a 
coherent Raman pulse technique using a laser system 
phase-locked on the two-photon hyperfine resonance 
\F = 1, m F = 0) -> \F = 2, m F = 0) (see Fig. 3a). Per- 
forming a 2tt pulse of duration 40 /is on one side of the 
condensate leaves the population effectively unchanged, 
but introduces a phase difference of ir compared to the 
unperturbed part of the condensate. In a small region 
around the edge of the mask I « £, a transfer of the 
population to the \F = 2, mp = 0) state occurs (Fig. 3b). 



By employing a step-like intensity pattern rather than a 
simple edge it is possible to vary the number of atoms 
transferred to the other hyperfme state by changing 
the width of the intermediate step. The mask pattern 
used in the experiment described here resulted in a 
bright component population Nb = 0.08A/" to t, where 
iVtot is the total number of atoms. A time series of the 
propagation of such a dark-bright soliton is shown in 
Fig. 3c. After a short time of flight of 9 ms the atoms 
are first exposed to a light pulse resonant only with 
\F = 2). After another 2 ms the \F = 1) atoms are 
subsequently imaged. The dynamics of the dark-bright 
soliton could be followed for more than 2s as seen in 
Fig. 3d. We observe an oscillation with a frequency of 
Q dh = 2tt x (0.90 ± 0.02) Hz = 0.24 x Q much smaller 
than the frequency of the corresponding dark soliton. 
An expression for the oscillation frequency can be given 
employing the equation of motion for very strongly 
populated dark-bright solitons [35] and leads to 
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where R z denotes the radius of the BEC along the axial 
direction and ot Zinax is an amplitude-dependent numerical 
factor which takes into account, that the total potential 
experienced by the bright component depends on the in- 
homogeneous density of the dark component. For the 
observed values of iVb, R z and N tot we calculate an oscil- 
lation frequency of 1.27 Hz which is in very good agree- 
ment with the observed value. 

Another spectacular feature that can be extracted from 
this measurement is the interaction of a dark soliton 
with the much slower dark-bright soliton. Owing to the 
method of initial state preparation, an extra dark soliton 
is always generated in addition to the dark-bright soliton. 
As shown in Fig. 4 the dark soliton propagates in the op- 
posite direction as compared to the dark-bright one and 
oscillates back with the same frequency as a dark soli- 
ton in an unperturbed experiment. After 120 ms it thus 
approaches the position of the dark-bright soliton which 
has only moved very little due to its much smaller os- 
cillation frequency. The dark soliton is reflected off the 
dark-bright one comparable to a hard- wall reflection and 
moves back. To our knowledge this is the first observation 
of collisions of different types of matter wave solitons. In 
conclusion we have realized long-lived solitons in 87 Rb 
BEC and observed soliton oscillations for the first time. 
Via a combination of a local Raman transfer and phase 
imprinting method we could realize two-component soli- 
tary excitations called dark-bright solitons which clearly 
exhibit very slow oscillatory dynamics in a trapped BEC. 
As a first striking example for soliton interaction the re- 
flection of a dark soliton bouncing off a filled soliton could 
be observed. These experiments pave the way for fur- 
ther studies on solitons and soliton dynamics in ultracold 
quantum gases. 



Methods 

Creation of solitons 

We create Bose-Einstein condensates of 87 Rb atoms 
via trapping of up to 5 x 10 9 atoms in a 3D magneto- 
optical trap (MOT), sub-doppler cooling of these atoms 
and subsequently transferring them into a magnetic trap. 
We evaporatively cool the atoms slightly above the criti- 
cal temperature T c for Bose-Einstein condensation within 
20 s. Afterwards we superimpose a crossed dipole trap 
realized by a Nd:YAG laser beam focused to a waist of 
35 /im and a perpendicular Ti:Sa laser beam (830 nm) fo- 
cused to a waist of 125 /im. The atoms are loaded in this 
dipole trap and further cooled evaporatively by smoothly 
lowering the optical power of the dipole trap beams un- 
til we end up with an almost pure BEC of 5 — 10 x 10 4 
atoms. We reduce the dipole trap power as much as 
possible in order to ensure smallest temperatures of the 
BEC. The trapping frequencies for this configuration 
have been determined by several independent approaches 
and read 2i\ x (6.3, 85, 133) rad s _1 . At this stage the con- 
densate consists of atoms in the \F = 1, mp = — 1) state 
and can be transferred to any other state or superposi- 
tion of states via rf- or microwave-pulse or -sweep tech- 
niques. The state preparation is carried out at a suffi- 
ciently large magnetic offset field to avoid undesired spin 
mixing dynamics [36 . Moreover, the creation of dark- 
bright solitons demands a spatially selective transfer to 
the \F — 2, m^? = 0) state which is accomplished by the 
use of a phase-locked Raman laser system with a rela- 
tive phase error of not more than 0.44 rad. The optical 
transfer and phase imprinting is achieved by imaging a 
computer-generated pattern displayed on a spatial light 
modulator with a pixel size of 8 /mi onto the BEC with an 
optical resolution of better than 2 /im. This is achieved 
via a high-quality imaging optics with a magnification of 
1/10. This optics is also used for the detection of BEC 
then yielding a magnification of 10. After creation of 
BEC in the optical dipole trap and the application of the 
phase imprinting light pattern via the SLM for a dura- 
tion of 40 /is, we allow the BEC to evolve in the trap for 
a variable time r ev oi- Subsequently, we switch off the op- 
tical trapping potential within 1 /is. In the experiments 
presented here we take an absorption image of the ex- 
panded atomic cloud after a time-of- flight of 11 ms. Note 
that this is a destructive detection technique. For the 
measurements presented here, each data point therefore 
corresponds to a new realization of BEC, phase imprint- 
ing, evolution time and detection. 



Determination of soliton parameters 

We have determined soliton parameters such as posi- 
tion, width and amplitude from 2D-fits to the absorption 
images. The employed function consists of a Thomas- 
Fermi like density distribution for the BEC modulated 



with individual solitons basically given by the square of 
Eqn. [2] The bright component has been fitted employ- 
ing a function derived from the wavefunction of the bright 
component of the dark-bright soliton [35 
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FIG. 1: Figure 1 — Principle of dark soliton gener- 
ation, a Optical setup: A spatial light modulator (SLM) is 
used to imprint a phase step by exposing part of the con- 
densate to a far-detuned laser beam, b Theoretical curve of a 
dark soliton's density |^| 2 (^) and phase <j) ( ), as described 
by Eqn. 2. c A typical absorption image of the condensate, 
taken directly after preparation of the soliton and a subse- 
quent free expansion of 11ms. Optical density is color- and 
height-coded for better visibility, d Integrated column den- 
sity (•) of the data in c together with a fit to the data (~). 
e Image of a soliton after an evolution time of 2.8 s. 




H 1 

0.75 

c 0.5 

o 

S 0.25 



8. 
§-0.25 
I -0.5 

to 

-0.75 
-1 






I 1 

'0.75 

0.5 

0.25 



-0.25 

-0.5 

-0.75 

-1 



4 ■ 
ft. 






▲ 
It 4 


A 4 

• 










J 1< 

' 1 






A * 


\ 4 


: 




A ▲ 


\ A " 




y 








\ * 




/ 


•ft 


k ▲ 


▲ 
• 




A 
• 


i 


• 


i 


4 

▲ 
* 



50 100 150 200 250 
time [ms] 



250 500 750 1000 1250 1500 1750 
time [ms] 



2000 



FIG. 2: Figure 2 — Dark soliton oscillations in a 
trapped BEC. a A set of absorption images showing the 
soliton position at various times after phase imprinting. The 
soliton propagates to the right and is reflected off the edge 
of the condensate after t ~ 80 ms. The corresponding evo- 
lution time for each image is given in units of the oscillation 
period T. b Results of a numerical calculation solving the ID 
Gross-Pitaevski equation corresponding to our parameters in 
units of T are presented. Experimentally observed features 
like density modulations caused by a density wave on the left 
side of the condensate as well as the development of a tiny 
second soliton are reproduced, c Axial positions of the soli- 
ton (•) with respect to the center of mass and normalized 
to the width of the condensate. The oscillation frequency is 
Q, — 2tt x (3.8=b0.1) Hz. The positition of a second tiny soliton 
( ) as well as a sinusoidal fit (^) to the position of the soli- 
ton are shown. Each data point was obtained from a different 
experimental run. The scatter is due to small fluctuations in 
the preparation process. Errors in extracting the solitons po- 
sition from the individual images are typically less than 0.02 
and therefore not plotted. 



FIG. 3: Figure 3 — Creation and oscillation of a dark- 
bright soliton. a A local population transfer in the center 
of the trapped BEC is achieved by a coherent two-photon 
Raman process between the two hyperfine states F = 1 and 
F — 2 leading to the generation of a dark-bright soliton (b). 
c A set of double exposure absorption images showing the 
density distributions of the two components which undergo 
slow oscillations in the axial direction, d Time series of the 
axial positions of the dark (•) and bright (a) component of 
the soliton in addition to corresponding sinusoidal fits to the 
position. Note that the time scale is different by almost an 
order of magnitude as compared to Fig. 2c. For details of the 
first 175 ms see Fig. 4. 
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FIG. 4: Figure 4 — Collision of a dark and a dark- 
bright soliton. A detailed plot of the first 175 ms of Figure 
3 reveals the reflection of an extra dark soliton off the dark- 
bright soliton. The axial position of the extra dark soliton (■) 
is plotted together with a fit to the data (^). The mean e -2 
width of the bright soliton is indicated ( ). The reflection 
of the extra dark soliton is very close to that expected from 
a hard- wall reflection. The fit corresponds to a sine- function 
mirrored at the reflection time t r = 117 ms. 



